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Abstract 

The theory of nematic-smectic phase transition in the system of 
uniform semi-flexible chains with hard-core repulsion is presented. 
Both the general density-functional formalism the tube-model cal- 
culation show that the flexibility of the chains results in a strong 
first-order transition, in contrast to the common weak-cristallization 
scenario of the nematic-smectic transition in rigid rods. The calcu- 
lated spinodal volume fraction of the uniform nematic phase and the 
period of the modulation instability are consistent with recent exper- 
imental results. 

PACS numbers: 64.70.Md, 61.30.Cz, 61.41. +e 
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1 Introduction 



The nematic-smectic phase transition is among the most interesting phenom- 
ena in hquid crystals. This is one of a few examples of a partial breaking of 
translational symmetry: the system in smectic state has a layered structure, 
while the translational symmetry along the in-layer directions is preserved. 
It is believed that normally the physical origin of the transition to the smec- 
tic phase is the non-uniform architecture of the consistuent particles, such 
that the nematic-smectic transition can be interpreted as a microphase sep- 
aration of different parts of the molecules. Nevertheless, it is conceptually 
important to realize that a non-trivial structure of the constituent objects is 
not a necessary condition for the formation of smectic. Both analytical and 
numerical investigations show that an entropy-driven nematic-to-smectic 
phase transition occurs even in the system of uniform rigid rods 

Recent experiments with rod-like viruses have confirmed the predictions 
of such theories 0], They also showed that a dramatic change in behavior 
takes place due to the finite flexibility of viruses. In particular, the transition 
turns out to be of the strong first order rather than of weak-crystallization 
type typical for most smectics. In addition, the volume fraction at the tran- 
sition is considerably higher in the case of flexible chains than for rods. The 
chains are strongly localized within the layers, so that the periodicity of the 
smectic phase is nearly equal to the length of an individual stretched chain, 
while in the stiff-rod case it is longer than the length of a rod because of a 
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weaker localization. 

In a recent communication we have proposed a theory based on the 
tube-model for the description of the nematic-smectic phase transition in 
the system of uniform semi-flexible chains with hard-core repulsion. This 
model has been shown to capture all the experimentally-observed features of 
the transition. In this paper we elaborate our approach by relating the tube- 
model calculation to a general density-functional formalism. It is shown that 
the strong first order character of the transition and the equivalence of the 
smectic period to the chain length follow from a general form of the density 
functional of the system. 

The strategy of the work is as follows. In Section 2 we review the density 
functional approach to the nematic-smectic transition. In Section 3 we de- 
rive the general form of the free energy of semi-fiexible chains expressed as 
a functional of the density of their mid-points. The breakdown of the sta- 
bility condition (positivity of the inverse structure factor) yields the limit of 
stability of the nematic phase and the critical wave vector of the modulation 
bifurcation. In Section 3 we obtain the parameters of the density functional 
from the tube model of the nematic state. This allows one to determine 
the parameters of the transition for various chain lengths. The structure of 
the theory allows one to understand the difference between the behavior of 
semi-flexible chains and that of rigid rods. 
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2 Density— functional Approach to Nematic— 
Smectic Transition 



A powerful tool for the description of various types of crystallization, and in 
particular, of the nematic-smectic phase transition is the density-functional 
approach. Within such theories the free energy F of the system is parame- 
terized by the one-particle density, p(r). The thermodynamic stability (or 
at least meta-stability) of the spatially uniform (e.g., nematic) phase with 
respect to density modulation is controlled by the sign of the corresponding 
second variation of the free energy: 



Note that G^^ (q) has the meaning of the inverse structure factor in the 
nematic phase (G (q) = {6pq6p-ci))- The system becomes unstable with 
respect to the transition to a spatially modulated (smectic) state, when in- 



nematic-smectic transition can be successfully described within the weak 



sition, the local deviations of the density field are dominated by one or by 
several critical density waves, i.e. 




(1) 



equality (|I]) is violated at a certain finite wave vector, qg. In many cases the 



crystallization theory [|T0] . In this approach one assumes that near the tran- 



n 




(2) 



a=l 



4 



Here are the critical wave vectors. The spatial dependence of the n- 
component order parameter pa is supposed to be much slower than the critical 
density modulation itself. 

The Landau-type expansion of the free energy in terms of this order 
parameter is the essence of the weak crystallization approach: 

f./(l,p= + ^V + ^V + fEIVp.|^)d. (3) 

The translational symmetry of the free energy demands that any term in 
the above expansion is a combination of density waves with zero total wave 
vector. Note that a large third-order term in this expansion would result 
in a strong first order transition and thus would violate the requirement of 
smallness of the order parameter. Therefore, the weak crystallization theory 
is adequate only for the cases when the third-order term is zero or small. 
In particular, its applicability to the nematic-smectic phase transition is 
normally justified by the broken rotational symmetry of the nematic phase. 
Indeed, in the rotationally-symmetric case, the critical wave vectors would 
lie on the sphere of radius qo, and one could choose three critical density 
waves with zero total wave vector. Therefore, rotational symmetry implies 
that there is a non-zero cubic term in expansion (^. In the case of the 
nematic-smectic transition the degeneracy in the orientation of the critical 
wave vector is lifted, and one cannot construct a third-order combination of 
the critical density waves. Hence, the weak crystallization theory can usually 
be applied to this transition. 

As an example, consider the transition to the smectic state in the system 
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of perfectly aligned hard rods. The density-functional theory, developed for 
this system in reference 0, shows that its basic physics can be successfully 
described even in the second virial approximation: 



which yields the following simple expression for the inverse structure factor: 



Here we are interested only in the wave vectors parallel to the nematic direc- 
tor, e. The breakdown of the stability condition corresponds to the transition 
to the spatially-modulated smectic phase. Since the instability is dominated 
by a single density wave (with q ~ ±1.57re/L), the transition in the hard- 
rod system is of the weak-crystallization type. This means that the smectic 
modulation of hard rods is weak near the transition point. As a result, the 
normal-to-layer fluctuations of the rods are of the order of their length L, 
consistent with the fact that the corresponding period A* ~ 1.3L differs con- 
siderably from L. Although within the model of freely-rotating rods 
the phase transition turns out to be of the first order, it is so weak that the 
transition can hardly be distinguished experimentally from a second-order 
one. 
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2.1 Nematic— Smectic Transition of Semi— Flexible Chains 

Consider a system of worm-like uniform chains of length L with hard core 
diameter D. The single-chain Hamiltonian for a given conformation r(s) 
(0 < s < L) has the form 

^ ^ 

Here s is the coordinate along the chain contour and = {dr/ds) is a unit 
tangent vector. The parameter p is the persistence length of the chain in 
the isotropic phase. In the nematic phase, however, there are two different 
scales which play the role of persistence length [O]. One of them is the 



typical distance between "hairpins", the points at which the tangent vector 
changes its direction by bl80°. This length becomes exponentially large for 
high enough nematic order parameter and we assume here that it exceeds 
the chain length, i. e. the conformations of the chains are straight lines with 
only weak transverse fluctuations of the tangent vector about the nematic 
director. The correlations of these transverse fluctuations are determined by 
another length scale known as the deflection length, C,±- This scale, which is 
smaller then the bare persistence length p, determines the thermodynamics 
of the system (the free energy can be estimated as kT per chain segment of 
length 

We now have to express the free energy as a functional of the density of 
the chains' centers, 
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p(r)= J: 5(r(L/2)-r) (7) 

chains 

If the chains were infinite, the conformational free energy averaged over scales 
larger then ^± would be a local functional of the volume fraction. For finite 
chains there is also a translational entropy contribution and a correction due 
to finite density of chain ends: 

Tp{chains) /. 

— ^ = y dr {p (r) Inp (r) + (0 (r)) + p^nd (r) /(^"^^ (0 (r))} (8) 

Because of the low density of chain ends, Pend, their effect is accounted for in 
the above expression by a term linear in Pend, coupled to some local function 
of the volume fraction. Since the chains are strongly stretched along the 
nematic director, z-axis, and do not form "hairpins", one can relate the 
local density of the ends p^nd (-2) and the local volume fraction of chains (z) 
to the density field p {z) (we will not consider any fiuctuations of these fields 
in the plane normal to z): 

Pend (z) = p {z+L/2) + p {z-L/2) (9) 

</.(;.) = ^ / p{z+s)ds (10) 

Due to these non-trivial relationships between the three fields p{r), (f> (r) and 
Pend(r), the above local free energy functional is becomes non-local when 
expressed in terms of a single density field, p. Such non-local properties 
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of the density functional are necessary for the description of the nematic- 
smectic transition. The semi-flexible chains are unique in that sense that 
this functional has a well-defined general form dictated by the property of 
locality on the mesoscopic scales (below L and above ,^^). 

In order to study the stability of the nematic state consider the second 
variation of the free energy at the fixed average volume fraction $ ( and the 
corresponding density p = $/Vo, where Vq = nD'^L/A is the volume of a 
single chain): 



^p{chains) 



fi fiend) 

Scf> (r) Scf> (r) + 2 



Spend (r) S(f) (r) 



(11) 



Performing the Fourier transform of all the fields under consideration ( along 
the ^-axis, neglecting their variations in other directions) and expressing 
these fields in terms of the density deviations, 5p, we obtain the following 
diagonalized free energy: 



gpichains) _Vo^ 5py6p. 



kT 



1 + A ($) 



L (1 — cosgL) sinqL' 
A($) {qLf 



(12) 



Here the dimensionless parameter A and the length A are certain functions 
of the average volume fraction, but not of the chain length or wave vector. 
Note the universahty of the (g, L)-dependence of the inverse structure factor 
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for semi-flexible chains: 
(q) 



I ^2p{chains) 



kT 6p,6p-, $ 



A /^x / L (1 — cosgL) singL' 



(13) 

Before proceeding with the discussion of this general form of G^^, we present 
a simple model which yields the microscopic expressions for the parameters 
A (<l>) and A ($) appearing in ([T3|) . 



3 Tube— model Calculation 

The interactions of a chain with its neighbors can be modeled by confining 
it in an effective tube. If the system is dense enough, it can be viewed as a 
close-packed array of such tubes. This means that the average tube diameter 
in the vicinity of some point r, is D/ ^(f) (r), where (r) is the local volume 
fraction. Therefore, the allowed amplitude of fluctuations of a chain within 
the tube is A = D [\/ ^(j){r) — 1^. The corresponding conformational free 
energy of the chain can be evaluated in the Gaussian approximation, with 
the confinement imposed by a fictitious external field. The resulting value is 
^kT per correlation length = (2A)^^^p^/^ [|l^. Note that one can obtain 
this result, up to a numerical coefficient, from simple scaling arguments, since 
^ can be identified with a typical contour length between two refiections of 
the chain from the tube walls. We conclude that the conformational free 



10 



energy per unit volume is determined by the field (r) : 




3 kT 40 (r) 
16^x(0) vrD2 



(14) 



We assume that the dominant effect responsible for the smectic ordering 
is that the internal parts of the chains cannot occupy a "shadow" region in 
the vicinity of a free end. This is a particular realization of the correlation 



behind the edge of one rod can be filled only by the complementary end part 
of another one (see Figure la). However, in the system of semi-flexible chains 
the size of the "shadow" region can be reduced by appropriate readjustments 
of the conformations of the neighboring chains, as shown in Figure lb. The 
screening of the "shadow" region can be described by a screening length Is- 
The conformational free energy penalty for the creation of the free space near 
the edge of every chain is given by the product of the local transverse pressure 
n_|_ = (l)df^^°"'^/d(p — /(™") and the typical volume of the "shadow" region 
'jirlgD'^/A. Here 7 is a geometrical factor of order of unity. Approximating 
the shape of a typical "shadow" region with a cone, one obtains 7 = 1/3. The 
typical bending energy associated with the distortion of chain contour needed 
for the screening of the "shadow region" is kTD'^p/lg per chain involved. The 
screening length and the energy of the end defect is determined by the balance 
between the osmotic energy penalty and the bending energy, i.e. they can 
be obtained by minimization of the following free energy: 



hole effect, |T^. In the extreme case of perfectly aligned rigid rods the space 




(15) 
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Here Z is the effective number of tlie distorted chains (fortunately, the final 
result is nearly insensitive to the choise of this parameter). The minimal 
value of the above free energy is given by 

^ I2^nf (iTrfV* ^ 2^1 f 1/^ _ (16) 



2 

The correspondong optimal screening length is 

/, = 2(2D)^/^pi/3 ^ 21* (l/^)-!^ (17) 

Here 1* = {2D)'^^^ p^^^ is the fundamental length scale of the problem. Up 
to a ^-dependent factor, it determines both the screening length and the 
correlation length ^± discussed earlier. For practical purposes, the screening 
length can be taken roughly equal to /*. Note that our estimate of the energy 
of the end defect is somewhat different from the one described in ref. [0], 
which is based on elasticity theory. Elastic approach neglects the energy 
penalty associated with the non-zero osmotic pressure H^, which in our case 
turns out to be considerably stronger then next-order elastic corrections. 
Nevertheless, the very notion of the effective attractin between chain ends 



discussed in the present paper is conceptually close to the one in ref. [|T4|. 

Depending on the total chain length, one can distinguish between two 
qualitatively different limiting regimes. If L <^ /*, the screening effect is not 
significant and the chains can be considered as rigid rods. Here, we consider 
the opposite limit, L 3> /*, when the total volume fraction of the "shadow" 
regions is low and one can neglect their overlap in the spatially-uniform 
nematic phase. 
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Summarizing the above calculation one can write down the total free 
energy of the system: 

= f / dr (l^oP (.) In p (.) + 

i (l/V^-l)' 



5/4 



(18) 



-(i/V^-r 

The first term accounts for the translational entropy of the chains, the second 
one represents the bulk conformational free energy of infinitely long chains 
and the last term is due to the end anomalies. One can easily verify that 
our model free energy has the general form obtained in the previous section, 
(^. There is a one-to-one correspondence between the three discussed con- 
tributions to the free energy and the three terms in the following expression 
for the inverse structure factor in the uniform nematic phase, as a function 
of the average volume fraction $: 



G [q) = 



2L (1-cosgL) 

- + 7 



(19) 

This expression has exactly the same general structure as has been derived 
in the previous section. The uniform nematic is stable (or at least metastable) 
with respect to the transition to the spatially-modulated (smectic) state as 
long as the calculated inverse structure factor is positive. The end effect 
contribution is the only term in expression (|Ty|) which may be negative (due 



to the sign-changing factor sin (qL) / qL) . For most wave vectors, however, 
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this term can not change the overall sign of the structure factor because of 
the dominant positive conformational contribution, which contains the large 
factor L/lg. This is related to the fact that the "shadow" regions are screened 
in the many-chain system, and the corresponding end effect is just a small 
correction to the conformational free energy. This correction is important 
only in the vicinity of the zeros of the expression (1 — cos qL) / (qL)"^, which 
determines the g-dcpcndence of the bulk conformational contribution to the 
inverse structure factor. 

We conclude that the modulation instability in the system is expected 
only for nearly L-periodic density waves, which are the zero modes of the 
conformational term. One can expand the inverse structure factor in the 
vicinity of such wave vectors, 27m/L {n — ±1, ±2, ...): 



(Sq) 



Vo I. 5 1 



1 



5q=q-2-Kn/L $ 16 ^ /i / A^'^ 



2)1 a \2Tin J \27rn^ 



$(l/\/$-l)' 

Its minima are the candidates for the critical wave vector of the modulation 
instability: 

g: = ^(l-3V$7^), n = ±l,±2,... (20) 

The control parameter $ at which G^^ (g*) = turns out to be independent 
of n (up to a cut-off nmax ~ L/l* where the small-(5g expansion becomes 
inadequate) : 
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The existence of the family of critical wave vectors ( pOD which differ only 
by integer multiplier is the signature of the first-order phase transition to 
the smectic state. Indeed, unlike the case of a single dominating density 
wave typical for most smectics (e.g. for rigid rods), this degeneracy en- 
ables one to compose the third-order combinations of the critical modes 
{ip (q*) ijj (gj^) ijj (qf), n + m + I = 0), which contribute to the term {Sip (z))^ 
in the density-deviation expansion of the free energy, (|]). The non-zero cubic 
term in the Landau expansion is known to result in a first-order transition. 

The common period of all the critical modes is 

X* = L + 3V¥-fl'{^*) (22) 

In a general case this critical period can differ from that of the equilibrium 
modulated phase. Nevertheless, the above result suggests that the period 
of the smectic A nearly coincides with the chain-length L, and that the 
small correction A — L is of order of Is — I*. This correction determines 
both the typical gap between well-formed smectic layers and the typical 
longitudinal fluctuations of the chains in these layers. This is consistent 
with the observations reported in ref. |^, for which the calculated length 
/* ^ lOOnm is of the order of the measured correction to the period A — 
L ~ 50nm. The effective hard core diameter in the experiments can be 
estimated as the inter-chain separation at which the electrostatic repulsion 
becomes of the order of kT. This length depends on the ionic strength, 
so it was possible to change it independently of the particle density. Since 
I* = {2DY/^p^f^ depends on the hard core diameter D, the measured points 
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on the phase diagram (concentration-ionic strength) can be transformed to 
the coordinates of our theory($ - L/l*). The calculated spinodal volume 
fraction, Eq.(pTD is in an agreement with the experimental value, which is 
about 0.75 for L/l* in the range from 4 to 10. Note that the chain-length 
dependence of the critical volume fraction is rather weak, as shown in Figure 
2. 

4 Discussion and Conclusions 

We now compare the results obtained for semi-flexible chains with those 
for perfectly-aligned rigid rods. Consider the inverse structure factor of 
perfectly-aligned rigid rods, Eq. (|): 

The first term here is due to translational entropy and the second one is 
the excluded volume contribution, which is essentially the "shadow" re- 
gion end effect. Naturally, the rigid-rod structure factor does not contain 
the conformational contribution, which dominates the similar expression for 
semi-fiexible chains. Hence, the "shadow" region effect, which drives for the 
nematic-to-smectic transition is no longer a small correction. As a result, 
the transition in the system of rods takes place at a lower volume fraction, 
which is $* ~ 0.57 within the second virial approximation, Eq. Taking 
the higher virial terms into account changes this value to 0.36. If the rods 
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are freely rotating, the critical volume fraction depends on the length-to- 
diameter ratio. For long enough chains {L/D ^ 10) the transition volume 
fraction reachs a "universal" value $ ~ 0.46 Q]. In the case of consid- 
erable flexibility of the "molecules", the critical (spinodal) volume fraction 
given by expression (^) is not constant even for high length-to-diameter 
ratio, since the relavant parameter here is -L//*, rather than L/D (see figure 
2). 

The typical behavior of the inverse structure factors for rigid rods, Eq. 
(^, and semi-fiexible chains, Eq. (JT^, are depicted in figure 3. Unlike 
the case of semi-fiexible chains {L ^ /*), in which the g*-degeneracy of 
the bifurcation point results in a strong first-order transition, the deepest 
minimum of the rigid-rod inverse structure factor determines a single critical 
wave vector of the modulation instability. This explains why the nematic- 
smectic transition for rigid rods is much softer then for semi-fiexible chains. 

Another important implication of the theory is the effect of polydisper- 
sity. Since the inverse structure factor describes the effective two-body in- 
teractions, expressions (|TP|) and can be extended to the polydisperse case 
by replacing a single parameter L with the mean length of two interacting 
particles (Li + L2) /2, and averaging over the distribution of lengths. The 
polydispersity acts against the modulation instabihty, because it reduces the 
depth of minima of the inverse structure factor. Hence, the critical volume 
fraction is expected to increase with polydispersity for both chains and rods. 
However, while the smectic phase can form in systems of rods with quite 
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broad distribution of lengths [|I5| , the typical deviation of the chain-length, 
which completely suppresses the transition in systems of chains is of the or- 
der of /*. Thus, in order to observe the smectic phase in the system of chains 
the distribution of lengths has to be very narrow. 

It should be noted that although we have shown that the above behavior 
follows from a very general form of the density functional for semi-flexible 
chains, our formalism cannot be directly applied to the case of long-range 
(hexagonal) in-layer structure. In that chain end becomes a real 

topological defect, and there is no reason to expect that the corresponding 
energy penalty is finite, i. e., that the corresponding contribution to the total 
free energy, Eq. (^, is linear in the chain-end density. 

In summary, the theory of nematic-smectic phase transition for uniform 
semi-flexible chains with hard-core repulsion has been developed. Similarly 
to the case of rigid rods, the transition is driven by the "shadow" region end 
effect. The difference is that due to the finite flexibility of the chains this 
effect is screened and the size of the empty space near a free end is limited by 
the screening length ~ Z* -C L. The presence of "shadow" regions yields 
just a small correction to the conformational free energy which stabilizes the 
spatially uniform nematic state. As a result, the spinodal volume fraction 
for flexible chains is much higher than for rigid rods. This trend is consistent 
with the experiments as well as with the recent theory of weakly flexible 
rods HTG). The modulation instability of the nematic state can appear only 
for nearly L-periodic density waves, which are the soft modes of the bulk 
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conformational free energy. Therefore, unlike the rigid-rod case, the period 
of the smectic phase near the transition point almost coincides with the 
chain length L. Another important difference is that in the case of semi- 
fiexible chains the point of modulation instability is highly degenerate in 
the critical wave vector resulting in a strong first-order phase transition. 
The theory also implies that the flexibihty of the chains results in higher 
sensitivity of the transition to polydispersity. The agreement of the theory 
with existing experimental data confirms that it captures the basic physics 
of the phenomenon. 
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Figure Captions. 



Figure 1. The "shadow" region (dashed) in the case of rigid rods (a) and 
semi-flexible chains (b). Note the screening oi the "shadow" by neigh- 
boring chains in the latter case. 

Figure 2. Spinodal volume fraction of uniform nematic $* as a function 
of reduced chain length L/l*. Solid line corresponds to the theoretical 
result obtained for semi-flexible chains {L/l* ^ 1), and the dashed one 
is an interpolation of the crossover to the rigid-rod limit {L/l* <^ 1). 
The geometrical factor 7 is taken to be 1/3, as is suggested in the text. 
Experimental points (diamonds) are taken from ref. [5]. 

Figure 3. Typical inverse structure factors (in arb. units) of the nematic 
phase for rigid rods (dashed line ) and semi-flexible chains (solid line) . 
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